Room- Temperature Superfluidity in Graphene Bilayers 
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Because graphene is an atomically two-dimensional gapless semiconductor with nearly identi- 
cal conduction and valence bands, graphene-based bilayers are attractive candidates for high- 
temperature electron-hole pair condensation. We present estimates which suggest that the 
Kosterlitz-Thouless temperatures of these two-dimensional counterflow superfluids can approach 
room temperature. 
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Introduction — Electron-hole pair (exciton) condensates 
were first proposed [Ij [2] as possible ordered states of 
solids more than forty years ago but have proved difficult 
to realize experimentally. Progress has been made re- 
cently with the discovery [21 13] of equilibrium exciton con- 
densation below T IK in the quantum Hall regime, the 
identification's] of spontaneous coherence effects in cold 
optically excited exciton gases, and studies of dynamic 
condensation fB] of polaritons in non-resonantly pumped 
optical microcavities. In the weak-coupling limit exciton 
condensation is a consequence of the Cooper instability [2] 
of solids with occupied conduction band states and empty 
valence band states inside identical Fermi surfaces. Bi- 
layer exciton condensates are counterflow superfluids 
with unusual electrical properties |iH71 [51 [TU[ [TT] which 
have so far been studied experimentally mainly in the 
quantum Hall regime. In this Letter we point out that 
superfluidity is likely to persist to remarkably high tem- 
peratures in graphene based bilayers. Graphene is a par- 
ticularly attractive candidate for room temperature bi- 
layer exciton condensation because it is atomically two- 
dimensional, because it is a gapless semiconductor, and 
because its two-dimensional massless Dirac band struc- 
ture implies nearly perfect particle-hole symmetry and 
stiff phase order. 

We consider a system with two graphene layers em- 
bedded in a dielectric media and gated above and below 
as illustrated schematically in Fig. [l] Each layer has two 
Dirac-cone bands centered at inequivalent points in its 
Brillouin-zone. The top and bottom gates can be used to 
control the electric fields E^r^t both above and below the 
bilayer. When the two fields are equal the bilayer is neu- 
tral, but charge is transferred from one layer to the other. 
The Fermi level lies in the graphene conduction band of 
one layer (the n-type layer) and in the valence band of 
the other layer (the p-type layer) . The particle-hole sym- 
metry of the Dirac equation ensures perfect nesting [12 
between the electron Fermi spheres in the n-type layer 
and its hole counterparts in the opposite layer, thereby 
driving the Cooper instability. The condensed state es- 
tablishes spontaneous long-range coherence between the 
two graphene layers. 




FIG. 1: (Color online) Left: Schematic illustration of a 
graphene bilayer exciton condensate channel in which two 
single-layer graphene sheets are separated by a dielectric 
(Si02 in this illustration) barrier. We predict that electron 
and hole carriers induced by external gates will form a high- 
temperature exciton condensate. Right: The two band model 
in which the two remote bands indicated by dashed lines are 
neglected. 

Our main interest here is in providing an esti- 
mate of the maximum possible Kosterlitz-Thouless (KT) 
temperature Tkt of these two-dimensional counterfiow 
superfiuids^. We use a two band model [13^ in which the 
occupied valence band of the n-type layer and the empty 
conduction band of the p-type layer are neglected. Our 
Tkt estimate is constructed from mean-field (Hartree- 
Fock) theory calculations [T3j of the temperature depen- 
dent phase stiffness of the ordered state. 

Our main result is the normal to superfluid phase 
boundary depicted in Fig. [2| The KT temperature is 
plotted as a function of the separation between the lay- 
ers d and the electric field E^^t outside the bilayer. We 
estimate that superfluidity can survive at room temper- 
ature under favorable experimental conditions. The non- 
monotonic dependence of Tkt on d at fixed Eg^t follows 
from a competition between the increasing carrier den- 
sity and the decreasing strength of interlayer electron- 
hole interactions with increasing d. At small d the phase 
stiffness is limited by the carrier concentration, which in- 
creases with d. At large d, the KT temperature is limited 
by the same fermion-entropy effects which are responsible 
for the Bardeen-Cooper-Schrieffer (BCS) critical temper- 
ature of weak-coupling superconductors. 
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FIG. 2: (Color online) Normal to superfluid phase diagram 
showing the dependence of the critical temperature Tc in 
Kelvin on the distance between layers d in nm and external 
bias electric field E^xt in V/nm. 



where e is the dielectric constant characterizing the em- 
bedding media, and 4>k = tai,n~^ (ky/kx)- The corre- 
sponding matrix element of the interlayer interaction is 
— V'^^pexp(~|fc — p\d). All energies are measured 



relative to the Dirac-point chemical potential of the bal- 
anced bilayer[TB]. Note that each spin and valley pairs 
independently and that electron-hole condensation is in- 
different to spin- valley space rotation in either layer. 

The interaction strength in a graphene monolayer is 
usually characterized by the dimensionless effective fine 
structure constant, a — e'^/ehv. This constant natu- 
rally appears in Eq. ([2| if energies and momentum are 
expressed in units of hvk^' and fep"' respectively. Here 
hvkp'^ = Vg/2 is the band Hamiltonian Fermi momen- 
tum. The strength of the interlayer interaction is deter- 
mined by a and by kp^d. 



Two-band mean-field theory — In the band eigen- 
state representation our band-Hamiltonian is Ti^ = 

-T,k,a',a'i,a'^kT^'a(^k,a, whcre Efc = Vg/2 - kvk, V 

is the band quasiparticle velocity, t is a Pauli matrix 
vector which acts on the which layer pseudospin, and 
Vg = eEextd is the gate induced potential difference be- 
tween the two layers. 

Spontaneous interlayer coherence is induced by inter- 
layer Coulomb interactions. In the mean- field descrip- 
tion the interlayer interaction reorganizes the low-energy 
fermionic degrees of freedom into quasiparticles which 
are phase coherent linear-combinations of the single-layer 
states. The mean- field theory Hamiltonian can be writ- 
ten in the following form fl5J: 



(1) 



Because of the model's particle-hole symmetry A° van- 
ishes. The pseudospin effective field Afc in Eq. ([ij solves 
the following self-consistent equation: 
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where A is the area of a graphene layer, Aj^ = (A^, A^), 
n is a unit vector parallel to A^, g = 4 accounts for 
the spin and valley degeneracy, and fd{x) = tanh(a;/2T) 
is the difference between the occupation numbers of the 
negative energy and positive energy quasiparticles. The 
Coulomb matrix element of the intralaycr interactions in 
the eigenstate basis is 
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FIG. 3: (Color online) Mean-field theory energy bands for 



a = 1, T = and 4"'^ = 1 



= Vg/2.) Note 



that kp < k„ because E^xt is screened. 



Interestingly, the self-consistent equations (|2| admit 
solutions with non-zero chirality J of the gap function 
A-*-: A^ = A^ (cos(J0fc),sin( J0fc)). However, the crit- 
ical temperature of a state with non-zero chirality is 
higher than that of the corresponding Tc of the zero chi- 
rality ground state so these solutions are unlikely to be 
physically relevant. We focus on the J = solutions 
hereafter. 

In the normal state, there is no interlayer coherence 
so A-'- vanishes. The intralayer Hartree-Fock potential 
then follows from self-consistent solution for A^. The 
main effects of electron-electron interactions in this case 
are to increase the bare quasi-particle velocity [T7j and to 
screen the external bias voltage. Screening reduces the 
amount of charge transfer and therefore reduces the nor- 
mal state Fermi momentum. As illustrated in Fig. |3] the 
energy bands change qualitatively in the condensed state 
because interlayer interactions induce coherence between 
the two layers and open an energy gap. 
Linearized gap equation — The mean-field theory phase 
boundary between the normal phase and superfluid phase 
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is obtained by solving the linearized gap equation 

i^(Afc) = i^A4,pn^(Ap) (4) 
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obtained by linearizing Eq. ^ with respect to A^. The 
kernel 
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of the linearized gap equation is obtained by solving the 
self-consistent equation for A^ in the normal phase. The 
normal phase is stable provided that all the eigenvalues 
of M are smaller than one. By numerically evaluating M 
for various interlayer distances and external fields wc find 
the mean-field phase diagram T"^{d, Eext) (not shown). 
Phase stiffness — In two-dimensional superfluids the crit- 
ical temperature is often substantially overestimated by 
mean-field theory and is ultimately limited by entrop- 
ically driven vortex and antivortex proliferation at the 
KT temperature 
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We estimate Tkt by using mean-field theory to calcu- 
late the phase stiffness (superfluid density) Ps{T). In 
parabolic band systems, this procedure yields reasonable 
estimates of Tkt in both BCS and BEC limits. 

The phase stiffness is most easily calculated by eval- 
uating the counterflow current jq — {e/h)psQ at small 
exciton momentum Q. Put formally, wc evaluate the 
expectation value of the counterflow current operator 



(7) 



with the density matrix defined by the mean-field Hamil- 
tonian 



^MP =n- + Y^ (^iQ^i+Q/2,^Ck-Q/2,l + h.C.) , (8) 



where A^q is the finite momenta pairing potential. 

Placing Q along x wc find that A^ ^hQv cos (j)k and 
that 



hvk^^^^-^-UiAk)h. (Afc) 



(9) 



(A^ = Efc). This expression for jg has an ultravio- 
let divergence and fails to vanish in the normal state 
(A-'- 0). Both properties are pathologies of the Dirac 
model. When the two Fermi circles are shifted in op- 
posite directions at finite Q they are asymmetric with 
respect to the momentum-space origin. As a conse- 
quence an ultraviolet cutoff at some momentum mag- 
nitude yields a finite counterflow current. This current 



would vanish if the same calculation was performed us- 
ing a microscopic model with integrations over the full 
Brillouin-zone. As long as A-'^ is small compared to 
graphene's 7r-band width, a condition that is very eas- 
ily satisfied, the anomalous ultraviolet contribution to 
Ps{T) is identical in the normal and in the condensed 
states. It follows that the physical counterfiow current is 
related to the Dirac model counterflow current (j°) by 
jg = iqiA) — jg(A-'- = 0). Following this prescription, 
we conclude that the last term in Eq. ([9| can be neglected 
and flnd that 
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Note that the zero temperature phase stiffness. 
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is purely a normal state property just as in BCS theory. 
Indeed an identical result is obtained in the BCS theory 
of a parabolic band system when ps is expressed in terms 
of the Fermi energy. 

An alternative approach for estimating Ps(T) which 
also accounts for the intralayer interactions is to evalu- 
ate the density matrix in Eq. ^ using the self-consistent 
mean-field equations with finite pairing momentum. As 
explained above, the physical counterflow current is ob- 
tained by subtracting jg(A-'- = 0) from Jq. The 
KT temperatures which follow from this procedure and 
Eq. (|6| are depicted in Fig. |2] Since Ps{T) is a de- 
creasing function of d it follows from Eqs.( 6|11 ) that 



In our calculations we find that this in- 



equality approaches an equality when k^d is small. Con- 
sequently, the increase in Tkt with d at small d in Fig. [2] 
simply follows the increase in E-p eEf,xtd/2. 
Discussion — The high transition temperatures we pre- 
dict deserve comment. They are larger than those of 
typical superconductors because condensation is driven 
by Coulomb interactions over the full band width, rather 
than by phonon-mediated interactions between quasipar- 
ticles in a narrow shell around the Fermi surface. In 
this sense exciton condensation is more akin to ferro- 
magnetism, which is also driven by Coulomb interactions 
and can survive to very high temperatures. The temper- 
atures at which exciton condensation can be achieved 
in graphene bilayers are immensely higher than those 
which might be possible in semiconductor bilayers be- 
cause more carriers can be induced by external electric 
fields when the semiconductor has no gap, because the 
Fermi energy increases more rapidly with carrier density 
for Dirac bands than for parabolic bands, and because 
graphene layers are atomically thin - eliminating the layer 
thickness effects which substantially weaken Coulomb in- 
teraction in semiconductor quantum well bilayers. The 
numerical estimates reported in Fig. [2] were obtained us- 
ing a coupling constant appropriate for a Si02 dielectric. 
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The optimal dielectric for high exciton condensation tem- 
peratures should have a high dielectric breakdown field 
and a low dielectric constant, suggesting that a suitable 
wide-gap material is likely the optimal choice. 

Screening and other beyond-mean-field induced- 
interaction effects are difficult to describe. In the case of 
weakly interacting atomic gases induced interaction ef- 
fects can [18] either increase or decrease T^, depending on 
the number of fermion flavors g. For the present Coulomb 
interaction case, a static Thomas-Fermi screening ap- 
proximation with normal state screening wavevectors re- 
duces interaction strengths very substantially when spin 
and valley degeneracies {g = 4) are included. Mean- 
field-theory critical temperatures are reduced by a factor 
of ~ at small d in this approximation and by a larger 
factor at large d. On the other hand, when the screening 
wavevectors are evaluated in the condensed state there 
is little influence on Tkt at small kpd both because the 
large gap weakens screening and because Tkt is propor- 
tional to the Fermi energy and not to the interaction 
strength in this limit. All this leads us to suspect that at 
low-temperatures there is a first-order phase transition 
as a function of layer separation d between condensed 
and electron-hole plasma states, similar to the transitions 
studied experimentally |T5] in quantum Hall exciton con- 
densates and theoretically pUj in parabolic band bilayers. 

Because of spin and valley degrees of freedom, the exci- 
ton pairing we describe in this work is SU(4) symmetric; 
crudely speaking the system has four identical superflu- 
ids simultaneously. We therefore anticipate interesting 
consequences of slightly unequal electron and hole densi- 
ties, similar to anticipated effects associated with the spin 
degree-of-freedom in normal exciton condensates |21l . 
Because of this sensitivity, front and back gates which 
can control the electric fields on opposite sides of the bi- 
layer independently are highly desirable in experimental 
searches for graphene bilayer exciton condensation. 

Our finding that Tkt ~ Q.lEp in the limit of strong 
interactions between conduction band electrons and va- 
lence band holes is partially supported by experimen- 
tal studies [23] of fermionic cold atoms in the strong- 
interaction unitary limit. It implies that Tkt should 
approach room temperature when Ep is larger than 
^ 0.3eV {n larger than ~ lO^'^cm"^) and d is smaller 
than ~ 2nm. Experimental detection of spontaneous 
coherence through one of its characteristic transport 
anomalies |4j will be necessary to construct a quantita- 
tively reliable phase diagram. 
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